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ABSTRACT 

The distribution of , the number of OB stars per association or cluster, appears to follow a universal 
power-law form N~ 2 in the local Universe. We evaluate the distribution of V* in the Small Magellanic 
Cloud using recent broadband optical and space-ultraviolet data, with special attention to the lowest 
values of iV*. We find that the power-law distribution in N* continues smoothly down to JV,= 1. 
This strongly suggests that the formation of field massive stars is a continuous process with those in 
associations, and that the field stars do not originate from a different star formation mode. Our results 
are consistent with the model that field massive stars represent the most massive members in groups 
of smaller stars, as expected if the clustering law applies to much lower masses as is expected from 
the stellar initial mass function (IMF). These results are consistent with the simultaneous existence 
of a universal IMF and a universal clustering law. Jointly, these laws imply that the fraction of field 
OB stars typically ranges from about 35% to 7% for most astrophysical situations, with an inverse 
logarithmic dependence on the most populous cluster, and hence, on galaxy size and/or star formation 
rate. There are important consequences for global feedback effects in galaxies: field stars should therefore 
contribute proportionately to the volume of the warm ionized medium, and equal relative contributions 
by superbubbles of all sizes to the interstellar porosity are expected. 

Subject headings: stars: early-type — stars: formation — stars: statistics — galaxies: star clusters — 
galaxies: stellar content — galaxies: individual (SMC) 



1. INTRODUCTION 

It is commonly held that most massive, OB stars are 
found in stellar clusters, or associations, since their short 
(~ 10 Myr) lifetimes are not long enough to permit spa- 
tial dispersion from their natal companions. However, 
apparently-isolated, massive field stars are well-known to 
exist, including a class of runaway OB stars with unusually 
large (;> 30 kms -1 ) velocities. While runaway OB stars 
are generally believed to be kincmatically ejected from 
a parent OB association, the ordinary non-runaway field 
OB stars have been suggested to originate in a different 
mode of star formation from their counterparts in associa- 
tions. Suggested differences between clusters and the field 
in the stellar initial mass function (IMF) for high-mass 
stars (Massey 2002; Kroupa & Weidner 2003) support this 
possibility. This contribution explores the relationship be- 
tween OB associations and massive field stars. 

In recent years, it has emerged that the number of stars 
N* per cluster appears to follow a universal power law 
distribution: 



N(N*) dN* cx N~ 2 dN* 



(1) 



This has been found empirically for young, massive clus- 
ters (e.g., Hunter et al. 2003; Zhang & Fall 1999), super 
star clusters (Meurer et al. 1995), globular clusters (e.g., 
Harris & Pudritz 1994) and H II regions (Oey & Clarke 
1998). The apparent universality of this relation is emerg- 
ing as fundamental (e.g., Oey & Muhoz-Tuhon 2003; Oey 
& Clarke 1998; Elmegreen & Efremov 1997), similar to 
the constant power-law relation for the IMF. If individual, 



field OB stars have a fundamentally different origin from 
clustered OB stars, then this is likely to be manifest in 
the distribution of iV», near N* — 1, where iV* specifically 
counts massive stars only (here, m ;> 10 M Q ). In what 
follows, we examine the form of N(N*), for small AT», a 
regime that has not been investigated to date, to illumi- 
nate the relationship between massive star clustering and 
massive field stars. 

2. THE SMC SAMPLE OF OB STARS 

A study of OB star clustering properties and field stars 
requires high spatial resolution and essentially complete 
detection of the OB stars over a large area of a given 
galaxy. The Magellanic Clouds are optimal target galax- 
ies by virtue of their proximity and high Galactic latitude. 
Indeed, some of the earliest work on OB associations and 
massive star censuses was carried out for the Magellanic 
Clouds by Feast et al. (1960), Sanduleak (1969), Lucke 
& Hodge (1970), Azzopardi & Vigneau (1982), and Hodge 
(1985). More recently, the UBVR survey data of Massey 
(2002) for the Small Magellanic Cloud (SMC), together 
with UV photometry (Parker et al. 1998) from the Ul- 
traviolet Imaging Telescope (UIT), provide unprecedented 
depth, resolution, completeness, and broadband coverage 
for the massive star population over a large fraction of the 
SMC. The area observed by UIT is smaller than the optical 
survey, but covers the SMC bar, which includes most of the 
active star formation in that galaxy. Figure 1 shows the 
regions covered by the two surveys. SMC OB associations 
have previously been identified and catalogued by Hodge 
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Field and cluster OB stars 




Fig. 1. — Grayscale i?-band image of the SMC (Bothun & Thompson 1988) showing the area coverage for the UBVR survey (Massey et 
al. 2002; rectangular regions) and the UIT survey (Parker et al. 1998; circular regions). The optical fields arc 79' square, and the UIT fields 
arc 37' in diameter. North is up, east to the left. 



(1985) and Battinelli (1991), the latter from the photo- 
graphic OB census of Azzopardi & Vigneau (1982). Here, 
we re-examine the clustering properties of these massive 
stars in the SMC, using the modern datasets to first sys- 
tematically identify OB stars, and then to systematically 
identify groups or associations of them. We will especially 
focus on the statistical properties of the low- TV* regime. 

Parker et al. (1998) provide UIT photometry in the B5 
filter (A cff = 1615 A, AA = 225 A), which recently has 
been recalibrated by Parker, Cornett, & Stecher (2004, in 
preparation). The UBVR photometry of Massey (2002) 
was observed through a Harris filter set and transformed to 
the Johnson-Cousins system (Landolt 1992). Reddening- 
free Q indices were computed using the extinction law of 
Cardelli, Clayton, & Mathis (1989) with a total to selec- 
tive extinction ratio Ry = 3.1. Different combinations of 
filters were used to best select two different stellar sam- 
ples: Stars with initial masses ^ 10 M Q , corresponding to 
spectral types of B0 V, B0. 5 I and earlier, were selected on 
the basis of their optical photometry only; and stars with 
initial masses ;> 20 M Q , corresponding to spectral types 
of 09 V, B0 I and earlier were selected on the basis of both 
their UBVR and UIT photometry. We refer to the former 
as the OB sample, and the latter as the O-star sample, al- 
though we emphasize that neither sample consists of spec- 
troscopically confirmed stars. Our criteria are as follows, 
using an SMC distance modulus DM= 18.88 (Dolphin et 
al. 2001) and global average extinction A B = 0.53: 

m B < —4.2 + DM + A B 



< 15.21 
Qu-r.b-r < -0.84, 
Qb5-u,b-r < -1-45, 
Qb5-v,u-v < 0.03, 



(2) 
(3) 
(4) 
(5) 



where 
Qu-r.b-r 



(mu - m R ) 



Au - A/ 



(m B - m R ) 



Ab-A r 

= (mu - m R ) - 1.396 (m B - m R ) , (6) 

n ( \ Ab5 - Au 

Qb5-u,b-r = (m B5 - mu) -. t— (m B - m R ) 

A B — A B 

= (m B5 -mu)- 1-668 (m B - m R ) , (7) 

^1 I \ A B 5 - Ay 

Qb5-v,u-v = (m B5 - my) — (mu ~ my) 

Au — Ay 

= (m,B5 — my) — 2.718 (mu — my) . (8) 

The OB sample is selected from equations 2 and 3 only, 
while the O-star sub-sample meets all of equations 2 to 5. 

Table 1 shows all the stars having both UIT B5 pho- 
tometry, UBVR photometry, and spectral classifications 
listed by Massey (2002), excluding stars with objective 
prism classifications and those identified as uncertain spec- 
tral types. The ID number and classification from Massey 
(2002) are shown in Columns 1 and 2. Column 3 shows 
whether the star met the O-star selection criteria above. 
Our attempts to select only stars earlier than 09 V and 
B0 I were successful for 27 out of 31 stars; of the 4 not 
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selected, one had nebular emission affecting the photome- 
try and two others had weak metal lines (designation W 
in spectral type). The false positives were 8 out 26 stars; 
the majority of the false positives were weak lined stars 
(W), emission line stars, and peculiar stars. 

3. IDENTIFICATION OF OB ASSOCIATIONS 

Having identified two samples of 1360 OB and 382 O 
stars, we then used the friends-of-friends algorithm de- 
scribed by Battinelli (1991) to identify the associations of 
these stars. Our parameter A* refers specifically to counts 
of OB stars from these samples, identified as described 
above. Battinelli's algorithm adopts as the clustering dis- 
tance d s between associated members the value of d s that 
maximizes the number of clusters for A* > 3. All stars 
within d s of another member star are defined to be within 
the same group. (We do not distinguish between "group," 
"association," and "cluster" in this work.) Figure 2 shows 
the number of clusters A as a function d s near the peak in 
A for the OB sample (Figure 2a) and O-star sample (Fig- 
ure 26). For the OB sample, the characteristic d s = 28 
pc (97") is straightforward to determine, but for the O- 
star sample, it is apparent that A is less sensitive to d s . 
We adopt d s = 34 pc (117") for the O-star sample, cor- 
responding to the mean value over the considered range; 
the larger value of d s as compared to the OB sample is 
consistent with the lower stellar density for O stars. Al- 
though the clustering properties are relatively insensitive 
to the choice of d s within a few pc, we caution that larger 
values of d s will result in flatter slopes for the clustering 
law A(A*) (e.g., equation 1), as discussed below in §4. 
Our values of d s are much smaller than Battinelli's value 
of 60 pc; this is consistent with a much greater complete- 
ness in our data, yielding a higher density of OB stars, and 
therefore smaller clustering distance. 

Table 2 presents the groups having at least 3 stars that 
are identified from our OB sample. Column 1 shows the 
group ID number, columns 2 and 3 give the group centroid 
position in decimal degrees (J2000.0), columns 4 and 5 
give the group diameter D in arcmin and pc, respectively, 
and column 6 lists the number of stars in the group. Ta- 
ble 3 presents the groups identified from the O star sam- 
ple in the same way. The group diameters are defined as 
D = \{Aa + AS) following Battinelli (1991), where Aa 
and AS represent the maximum difference between mem- 
bers in RA and Dec, respectively, in degrees of arc. We 
also list the individual member stars for each group in Ta- 
bles 4 and 5, which are fully available in the on-line edition. 
The first two columns of Tables 4 and 5 give the RA and 
Dec of each star in decimal degrees (J2000.0), columns 3 
- 10 give the U BVR magnitudes and uncertainties from 
Massey (2002), and the last two columns give the star ID 
from Massey (2002) and our newly-determined OB group 
ID. Table 5 also lists the UIT B5 magnitudes from Parker 
etal. (2004) in columns 11 and 12. Note that these group 
identifications represent the results from two separate runs 
of the group-finding algorithm; thus, if a star belongs to 
both samples, it may belong to groups in both Table 4 and 
Table 5. 

Figure 3 shows the location of the group centroids com- 
pared to those of associations identified by Hodge (1985; 
solid black circles) and Battinelli (1991; dashed black cir- 
cles) . Our OB and O-star groups are indicated by the solid 



blue and red circles, respectively. The circle sizes corre- 
spond to the mean diameters D from Tables 2 and 3. For 
the Hodge objects, the diameters are taken to be the mean 
of the dimensions in RA and Dec given by Hodge (1985). 
There is good general correspondence between the posi- 
tions of our groups and these earlier catalogs; certainly at 
least as good as the correspondence between the Hodge 
(1985) and Battinelli (1991) identifications. It is apparent 
that the net tendency from our smaller d s is for our associ- 
ations to be smaller than those in the earlier catalogs, and 
to break up some of those associations into smaller groups. 
Figures 4a, c and 46, d show the size distributions for, re- 
spectively, the OB and O-star samples. It is apparent that, 
aside from individual stars, the peak in the distribution is 
around 10-15 pc, compared to 50 pc for both the Battinelli 
(1991) and Hodge (1985) catalogs. Whereas Hodge (1985) 
reports that the mean diameter of his sample is similar to 
the mean for the LMC sample, our smaller characteristic 
value is consistent with the lower luminosities of the SMC 
H II regions and star-forming regions compared to those 
of the LMC (e.g., Kennicutt, Edgar, & Hodge 1989). 

4. FIELD VS CLUSTERS 

Figure 5 shows the distribution in A* for the OB groups 
(panel a) and O-star groups (panel 6). Two slopes are 
shown fitted to the power-law distribution, weighted by 
the square root of the bin value, \J AT (log A*): the solid 
line shows the fit for the entire distribution, yielding 
-2.51 ± 0.29 and -2.33 ± 0.42 for the OB and O-star 
samples, respectively; and the dotted line shows the fit 
omitting the first bin, which corresponds to single stars, 
resulting in fits of -2.27±0.38 and -2.11±0.55 for the OB 
and O-star samples, respectively. The fits omitting the sin- 
gle stars agree with the power-law slope of -2 (equation 1), 
found for the A* distribution in a variety of systems, as 
discussed in §1. We again emphasize, as mentioned above, 
that the fitted slope has a dependence on the clustering 
distance d s : for large d s , more stars are drawn into the as- 
sociations, causing the resultant slope of the clustering law 
to be flatter relative to clustering defined by a small d s . 
Since we seek a characteristic value for the slope of AT (AT*), 
it is important to ensure that d s is in turn characteristic 
of the sample. For the d s in the extreme range of 20 - 
40 pc (see Figure 26) , the O-star sample shows a variation 
in fitted slope, omitting the field stars, of —2.79 ± 0.84 to 
— 1.85 ±0.52, respectively. 

It is apparent that the single stars are slightly greater 
in number than the power law for the remainder of the A* 
distribution, particularly in the OB sample. How signifi- 
cant is this excess in the isolated, "field" stars? We note 
that the magnitude of the excess may be a lower limit, 
since the slope of the A* distribution is expected to flatten 
slightly at the smallest values of A*, owing to statistical 
effects, as follows. We assume, as suggested above, that 
the OB stars counted in the value of A* represent only the 
most massive stars for a population of star clusters that 
are in reality populated by a mass distribution described 
by a conventional IMF. Since we identified only OB stars, 
our clusters are selected with the criterion that each clus- 
ter contains a star of at least mass m cut . The smallest, 
"unsaturated" clusters have a lower probability of having 
their maximum stellar mass m max > m cut , so they will be 
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Table 1 

Comparison of O-star criteria 



Cat ID a 


Sr> Tvr>p b 


Td pti t i fi pd c 

\y uL n > 


Cat TD a 


Sr> Tvnp b 


TH pti t ifi pd c 

\y 01 u j 


16828 


05 111(f) 


y 


53382 


B0 IIWW 


y 


38024 


05 V 


y 


54456 


B0 IW 


y 


5869 


07 V 


y 


54958 


B0 IIIWW 


y 


16056 


07 III 


y 


9079 


B1.5 V 


n 


16885 


07 III 


y 


9488 


Bl I 


n 


17457 


07 III 


y 


18700 


Bl lb 


y 


17927 


07.5Iaf+ 


y 


37732 


Bl IIW 


n 


19650 


07.5 III 


y 


42740 


Bl III 


n 


22837 


07.5 V 


y 


43807 


Bl V 


n 


40341 


07 III 


y 


44784 


Bl III 


n 


43724 


07 V 


y 


45114 


Bl III 


y 


45521 


07 If 


y 


45438 


Bl Ve 


n 


46035 


07 III 


y 


45809 


Bl III 


n 


52170 


07 V 


y 


50609 


Bl III 


n 


27731 


08.5 V 


y 


51009 


B1.5 III 


n 


43197 


08.5 V 


y 


53225 


Bl V 


n 


43734 


08 V 


n 


26901 


B2 IIIW 


y 


49580 


08.5 V 


y 


32907 


B2 II 


n 


50825 


08 Vn 


y 


50826 


B2.5 III 


n 


53373 


08.5 V 


y 


54818 


B2 I 


n 


6406 


09+NEB 


n 


24929 


B3 I 


n 


11777 


09 II 


y 


43844 


B3 I 


n 


40610 


09 V 


y 


44828 


B3 I 


11 


45677 


09 III 


y 


47028 


B3 I 


n 


47540 


09 III 


y 


49606 


B3 I 


n 


3459 


B0 IWW 


n 


45722 


B8 Vp 


y 


13075 


B0.5 III 


y 


50475 


B8 lab 


n 


13831 


B0.5 V 


y 


53084 


B8-A0 I 


n 


15503 


B0.5 V 


n 


5391 


B pec 


y 


15742 


B0 IIWW 


y 


21801 


B pec 


y 


18614 


B0 IIW 


n 


24914 


B extr 


y 


20656 


B0 I 


y 


40851 


B pxtr 


y 


23702 


B0 IW 


y 


10915 


AO la 


n 


41648 


B0.5 III 


n 


43215 


AO I 


n 


43686 


B0.5 V 


n 


48732 


AO I 


n 


43758 


B0.2 V 


n 


52992 


Al I 


11 


51575 


B0 III 


y 









a Catalog number from Massey (2002) 

b Spectroscopic type compiled by Massey (2002) 

c Identified by O-star criteria equations 2-5 
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Fig. 2. — Number of groups with N* > 3 as a function of clustering distance d s for the OB sample (panel a) and O-star sample (panel b). 



-71 .5 




Fig. 3. — Centroids and sizes of associations identified by Hodge (1985; solid black), Battinelli (1991; dotted black), our OB sample (blue) 
and our O-star sample (red). Note that all blue circles have a dot at the center, thus overplotted identical circles can be identified by small 
dots in the center of red circles. (RA and Dec given in degrees of arc.) 
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Table 2 

Groups from the candidate OB sample 



OB grp 


RA 


Dec 


D(min) 


D(pc) 


N„ 


20 


11.3029 


-73.3820 


1.38 


23.93 


4 


24 


11.3639 


-73.2441 


3.39 


58.90 


5 


29 


11.5932 


-73.1044 


1.45 


25.21 


4 


30 


11.5857 


-73.4439 


1.69 


29.37 


3 


31 


11.5817 


-73.2120 


0.38 


6.64 


3 


38 


12.0545 


-73.1842 


12.91 


224.68 


70 


40 


11.7724 


-73.2863 


0.81 


14.05 


3 


41 


11.8078 


-73.3658 


1.97 


34.36 


4 


45 


12.0111 


-73.3695 


1.92 


33.33 


4 


47 


12.2249 


-73.4337 


6.93 


120.56 


26 


64 


12.2657 


-72.0873 


1.75 


30.42 


4 


73 


12.4821 


-73.3470 


3.95 


68.70 


16 


74 


12.4058 


-73.0955 


0.85 


14.84 


3 


75 


12.4251 


-72.8051 


1.02 


17.81 


3 


86 


12.5262 


-73.2635 


0.90 


15.73 


3 


91 


12.6156 


-72.8806 


1.16 


20.18 


4 


98 


12.6898 


-73.4584 


1.47 


25.53 


6 


99 


12.7466 


-72.6852 


3.78 


65.71 


9 


100 


12.6862 


-72.7800 


1.99 


34.61 


5 


101 


12.7490 


-73.2845 


3.42 


59.56 


10 


104 


12.7620 


-73.1756 


1.14 


19.90 


4 


107 


13.1672 


-72.7186 


12.48 


217.09 


55 


110 


12.8384 


-73.3445 


2.35 


40.93 


7 


112 


12.8514 


-72.9475 


1.46 


25.40 


3 


115 


12.8844 


-72.6976 


1.20 


20.94 


3 


118 


12.9530 


-72.8498 


3.02 


52.62 


7 


119 


12.8984 


-72.4770 


0.91 


15.80 


3 


120 


12.9375 


-73.2152 


2.21 


38.45 


8 


121 


12.9762 


-73.3923 


1.59 


27.72 


3 


123 


13.1031 


-72.1594 


5.99 


104.24 


15 


125 


13.0386 


-72.6643 


2.14 


37.23 


6 


131 


13.0357 


-72.3948 


1.05 


18.26 


3 


137 


13.1042 


-73.3233 


1.57 


27.33 


4 


142 


13.3029 


-73.1472 


3.52 


61.32 


8 


150 


13.2429 


-72.9064 


1.64 


28.61 


3 


152 


13.2756 


-73.1909 


1.06 


18.42 


3 


158 


13.4218 


-72.5561 


1.30 


22.63 


3 


160 


13.4898 


-72.2771 


3.20 


55.66 


6 


170 


13.5753 


-72.3525 


3.01 


52.33 


9 


174 


13.6399 


-72.4905 


1.70 


29.66 


4 


195 


13.9908 


-72.3949 


1.86 


32.43 


4 


197 


13.9913 


-72.2634 


1.06 


18.43 


3 




1 A D9Q1 


79 


n 7i 

U. 1 1 




Q 
O 


211 


14.0923 


-72.0363 


0.82 


14.26 


3 


228 


14.3369 


-72.6518 


1.37 


23.86 


5 


233 


14.4076 


-72.3812 


2.52 


43.93 


5 


235 


14.4236 


-72.2069 


1.26 


22.01 


3 


237 


14.5052 


-72.5974 


1.74 


30.22 


5 


241 


14.5732 


-72.6575 


2.34 


40.66 


8 


244 


14.7403 


-72.1810 


5.74 


99.90 


42 


247 


14.6652 


-72.7249 


2.85 


49.57 


5 


248 


14.5854 


-72.3032 


0.78 


13.66 


3 


255 


14.6828 


-72.2829 


0.98 


17.03 


3 


256 


14.7037 


-72.3284 


1.59 


27.67 


3 


263 


14.8208 


-72.0852 


1.10 


19.18 


3 


265 


14.8466 


-72.3857 


1.06 


18.47 


3 


270 


14.8983 


-72.7153 


1.40 


24.44 


3 


277 


15.0082 


-72.2228 


2.09 


36.34 


4 


279 


14.9884 


-72.1773 


1.54 


26.76 


3 


294 


15.1744 


-72.5081 


1.87 


32.59 


4 
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Table 2 
( continued) 



OB grp RA Dec D(min) D(pc) N t 



303 
ouo 


1 ^ 9701 


79 931 7 


2 02 


3^ 1 & 
00. 10 


4 


31 ^ 


1 ^ 3803 


79 1 378 
1 z. 10 / 


1 38 
1 .00 


93 Qt\ 

zo.yo 


Q 
O 


39^ 

ozo 


1 ^ A91 £ 


71 QAR1 

— / ±.y4:Dl 


1 98 
1 . Zo 


99 91 
zz.zl 





097 

OZr I 


15 5244 


79 90Q1 


2 15 


37 AO 

O < .4:U 


r; 



331 


15 5112 


79 ^Q9& 


U.OO 


1 ^ 30 


Q 
O 


333 

000 


1 o.oooz 


79 093Q 
/ z.uzoy 


3 A7 

O . 4: I 


fiO A1 

UU.4: J_ 


7 


337 


iD.oooy 


79 OSQQ 
— / z.uoyy 


& A^ 

0.4:0 


1 A 7 1 O 

14: I . 1U 




3ao 

04:U 


1 ^ R331 

ID.OOOl 


71 Sfifi^ 


9 A1 

Z.4:± 


A1 Sfi 





344 


1 ^ 791 ft 


71 Q9fi1 


3 90 

O. ZrU 


^ 77 
OO. / I 


4 


3A C » 

04:0 




—79 7(^7^ 


1 Q9 

i.yz 


33 33 
OO.OO 


4 


3A7 

04: / 


1 ^ 7^R9 
10- / OOz 


79 A9^9 
— / Z.4ZOZ 


9 3^ 
z.OO 


AO 89 


a 



3A& 

04:0 


k 71 fin 


79 A&R7 


9 0.7 
z.u / 


3^ Q& 

00. yo 


4 


3^1 
OOl 


10. 1 uoo 


79 fi9A0 


1 Rfi 
l.OU 


32 41 


Q 
O 


354 


15.9167 


—71 QAfiQ 
1 _l . y4:uy 


3.32 


57.77 


g 


3^7 
00 i 


lO.OUOO 


79 3&77 


1 A3 

1 .4:0 


9A S 1 ^ 

Z4:.00 


Q 
O 


ouo 


1 fi ri99fi 


—79 73^ 

|Z. 1 OOO 


9 00 
Z.UU 


3A SO 

04:.0U 


Q 
O 


366 
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Table 3 

Groups from the candidate O-star sample 



O grp 


RA 


Dec 


D(min) 


D(pc) 


N* 


8 


11.9117 


-73.1385 


4.87 


84.71 


11 


11 


11.8196 


-73.2295 


1.11 


19.35 


3 


13 


12.0083 


-73.3533 


2.64 


45.86 


4 


19 


12.2742 


-73.2744 


3.57 


62.11 


6 


24 


12.2523 


-73.1099 


1.41 


24.49 


3 


26 


12.4694 


-72.6139 


3.03 


52.78 


4 


38 


12.6156 


-72.8806 


1.16 


20.18 


4 


41 


13.0895 


-72.7004 


13.76 


239.45 


57 


45 


12.9077 


-72.7988 


5.06 


88.03 


12 


50 


12.9830 


-72.8530 


2.30 


40.07 


4 


51 


12.9693 


-72.5435 


1.60 


27.84 


4 


62 


13.4522 


-72.5525 


2.19 


38.11 


4 


67 


13.5963 


-72.3412 


5.02 


87.31 


12 


68 


13.6438 


-72.4994 


2.52 


43.85 


5 


76 


13.8855 


-72.4820 


1.74 


30.28 


3 


78 


13.9908 


-72.3949 


1.86 


32.43 


4 


81 


14.2089 


-72.4833 


5.75 


100.05 


11 


84 


14.0923 


-72.0363 


0.82 


14.26 


3 


97 


14.4076 


-72.3812 


2.52 


43.93 


5 


102 


14.5581 


-72.6512 


1.87 


32.61 


4 


104 


14.7329 


-72.1834 


5.90 


102.61 


26 


107 


14.6211 


-72.7179 


2.00 


34.71 


3 


109 


14.6682 


-72.2968 


2.38 


41.43 


5 


121 


15.0359 


-72.2274 


1.72 


29.92 


3 


124 


15.0641 


-72.5572 


1.42 


24.74 


3 


130 


15.2501 


-72.2311 


2.68 


46.59 


5 


143 


15.5244 


-72.2091 


2.15 


37.40 


5 


145 


15.5583 


-72.0413 


1.55 


26.96 


3 


146 


15.7173 


-72.1075 


5.36 


93.19 


17 



Table 4 

Group membership for the candidate OB sample & 



ra 


Dec 


U 


U err 


B 


B err 


V 


V err 


R 


R err 


ID 


OB grp 


10.1171 


-73.5425 


14.01 


0.05 


14.96 


0.04 


15.00 


0.03 


14.96 


0.07 


107 


1 


10.1832 


-73.4063 


14.40 


0.05 


15.18 


0.04 


15.12 


0.03 


14.86 


0.07 


298 


2 


10.3880 


-73.4256 


14.24 


0.06 


15.15 


0.06 


15.28 


0.04 


15.30 


0.11 


1037 


3 


10.5417 


-73.2323 


13.44 


0.03 


14.42 


0.03 


14.60 


0.02 


14.69 


0.05 


1600 


4 


10.5515 


-73.3866 


14.26 


0.05 


15.19 


0.04 


15.15 


0.03 


15.04 


0.08 


1631 


5 



a The complete version of this table is in the electronic edition of the Journal. The printed edition contains 
only a sample. 



Table 5 

Group membership for the candidate O-star sample" 1 



ra 


Dec 


U 


U err 


B 


B err 


V 


V err 


R 


R err 


B5 


B5 err 


ID 


O grp 


11.1440 


-73.1597 


13.78 


0.05 


14.70 


0.04 


14.82 


0.03 


14.74 


0.07 


11.92 


0.05 


4424 


1 


11.2379 


-73.0130 


13.19 


0.03 


14.30 


0.03 


14.36 


0.02 


14.29 


0.04 


11.16 


0.09 


4922 


2 


11.2438 


-73.0213 


14.05 


0.05 


14.94 


0.04 


14.91 


0.03 


14.86 


0.07 


12.19 


0.05 


4949 


2 


11.3259 


-73.2564 


12.45 


0.02 


13.36 


0.01 


13.31 


0.01 


13.14 


0.02 


10.25 


0.08 


5391 


3 


11.3867 


-73.0773 


13.16 


0.04 


14.15 


0.04 


14.30 


0.03 


14.28 


0.06 


10.85 


0.06 


5718 


4 



a The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 
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Fig. 4. — Logarithmic size distribution for groups identified from the candidate OB and O-star samples are shown in the top row (panels a 
and b, respectively). The linear distributions for the OB and and O-star samples are shown in the bottom row (panels c and d, respectively). 
The first bin in Figure 4c has a maximum of 391. 




Fig. 5. — The distribution in N t for groups defined from the candidate OB sample (panel a) and O-star sample (panel b). Fits to the 
power-law exponent (—f3) are shown for the entire distributions (upper values) and omitting the first bin of single stars (lower values). Note 
that the slope fitted in logarithmic space is — (3 + 1. 
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progressively missing from a sample of clusters selected in 
this way (see Oey & Clarke 2003 for a detailed discussion 
of this effect). 

Figure 6 shows Monte Carlo models that demonstrate 
this effect. TV* is drawn from a power-law distribution of 
slope -2 (equation 1), with the individual stellar masses 
drawn from a Salpeter (1955) IMF having 10 < m < 
100 M Q . Figure 6a shows the distribution in for 10,000 
clusters having at least one star with mass m max > 20 M Q , 
and Figure 66 is the same, but selected for m max > 40 M Q . 
We sec a flattening in the distributions for small iV* , and 
the flattening is stronger in Figure 66, since the probability 
of drawing a 40 M Q star is less than for a 20 M Q star. 

Figure 7 shows the stellar density of all catalogued stars 
in the UBVR survey as a function of radius around the iso- 
lated, field massive stars found in the OB (solid line) and 
O-star (dashed line) samples. The errors are computed as 
y/n/nr for the total n stars within radius r, omitting the 
central field massive stars. As seen in Figure 7, it is ap- 
parent that the stellar density increases at small r. This 
is consistent with our suggestion that most massive field 
stars represent the most massive component of groups of 
smaller stars, as is expected from the stellar IMF. 

Another effect that can offset the statistical flattening 
in the form of the distribution for small iV* , is evolu- 
tion. In the most extreme situation, we consider that all 
the OB associations in the SMC were formed together in 
a single burst of global star formation. Staveley-Smith et 
al. (1997) find that the H I shells identified in the SMC 
appear to have a narrow age distribution around 5 Myr. 
We described above that the unsaturated, lowest-iV* clus- 
ters have a lower average stellar mass (e.g., Oey & Clarke 
2003). Since stellar lifetimes are longer for lower-mass 
stars, we therefore expect that these unsaturated clusters 
will tend to last longer, on average, than the statistically 
fully-sampled, or saturated, objects. Thus, we may expect 
that the distribution in iV* steepens with time. 

We constructed Monte Carlo simulations of such an ag- 
ing burst of clusters. Figure 8 shows the models for a 



population of clusters at 0, 15, and 25 Myr after their 
simultaneous formation. We used the same IMF param- 
eters are before, and the stellar ages are from the grid of 
Charbonnel et al. (1993) for SMC metallicity. Each of the 
model TV* distributions is fitted with a power-law slope, 
shown in Figure 8. We see that, while it is difficult to 
discern the slope steepening over the entire sample after 
15 Myr, and even at 25 Myr, it is apparent that the ag- 
ing effect is most pronounced in the smallest N* bins, and 
that the 7V*=1 bin becomes the most disproportionately 
enhanced. For continuous creation of the clusters, the net 
observed slope will be intermediate between the forms in 
Figure 6 and Figure 8. 

The modest observed excess of single stars is also likely 
to be caused in part by a contribution to their popula- 
tion by runaway OB stars. If runaway OB stars originate 
in associations, however, then this contribution must be 
small, since runaways correspond to about 3% and 20% 
of field early B and early O stars, respectively (Blaauw 
1961). Runaways will contribute primarily to the single- 
star population, while binary runaways represent <j 19% 
of all runaways (Gies 1987). Binary runaways are likely to 
be tight pairs that remain unresolved in the 0.67 pc px _1 
survey resolution. 

It is also possible that spurious selection of candidate 
stars, owing to the coarse photometric criteria, are a fac- 
tor. We showed in §2 that there may be a net of about 15% 
of the O-star sample that consist of spuriously selected B 
stars; while in principle, these should be distributed pro- 
portionately between clusters and the field, it is possible 
that crowding effects in clusters favor field selection. If 
so, owing to the longer lifetimes of B stars, they could 
contribute disproportionately to the field star population. 
However, we emphasize that this requires a significant vari- 
ation in the spatial distribution of only the spurious can- 
didates, between field and clusters. 

The most important factor in enhancing the field star 
population, however, is probably the strong variation in 
star formation density across the SMC. The highest star 
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Fig. 7. — Mean stellar density (arcsec -2 ) within radius r (arcsec) from isolated massive stars in the candidate OB (solid line) and O-star 
(dashed line) samples (10" = 2.9 pc). The O-star sample shows a higher background density since that sample is limited to the SMC bar 
region (Figure 1). 



formation density occurs in the SMC bar, where, as seen in 
Figure 3, most of the associations are located. Individual 
massive field stars, on the other hand, have a much more 
uniform distribution, as seen in Figure 9. The uneven star 
formation density distribution in this galaxy therefore fa- 
vors the field stars and enhances the N^— 1 bin for the OB 
sample in Figure 5. The effect is much less pronounced in 
the O-star sample, which is limited to the SMC bar region. 

Considering these effects, it is apparent that there is 
no strong variation or change in character of the power- 
law distribution seen in Figure 5 for the smallest values of 
N*. The empirical samples (Figure 5) are much smaller 
than those in the models and therefore have significantly 
poorer statistics in logiV(log TV*). Thus, the modeled ef- 
fects will be more difficult to discern in the data. Since 
our results are largely consistent with a single intrinsic 
power-law form for the clustering law, this suggests that 
the massive field star population simply represents an ex- 
tension of the massive cluster population extending down 
to iV* = 1. There is no evidence that the majority of 
massive field stars originate from a mode of star forma- 
tion that is different from those in associations. However, 
further studies of massive field and cluster populations in 
other environments are necessary to confirm the generality 
of these results. 



5. DISCUSSION 

For a universal IMF, the constant slope of the dis- 
tribution extending to the field stars has profound conse- 
quences for their global feedback influence in galaxies. We 
can now quantify the assertion that most massive stars 
form in associations: for an distribution given by equa- 



tion 1, the total number of OB stars is, 

iV,,up 

AT*,tot oc N * ■ N * 2 > ( 9 ) 

JV» = 1 

where -/V* :Up refers to the upper-limit, maximum cluster 
of stars. Equation 9 corresponds to a divergent harmonic 
series which can be approximated for large iV*, up as: 

AT*,tot°5 lnAT* iUp + 7 , (10) 

where 7 ~ 0.5772 is the Euler-Mascheroni constant. Thus 
the fraction of JV*= 1 field stars is (In -/V* iUp + j)^ 1 of the 
total iV^tot- For our OB and O star samples, respectively, 
logA^up ~ 2.0 and 1.8, yielding 19% and 21% fractions 
for the field stars. 

Counting the actual stars in our OB and O star sam- 
ples, we find that 374 and 91 candidates, respectively, 
had no massive companions within the clustering radius. 
These correspond to 28% and 24%, reflecting the excess 
found above. These fractions are about a factor two lower 
than the finding by Parker et al. (2001) that over half 
of their i/IT-selected candidate O stars are outside cat- 
alogued association boundaries in the Large Magellanic 
Cloud (LMC). The results for both the SMC and LMC 
may be odds with the results of Massey et al. (2002), 
who found a much steeper IMF slope for field vs cluster 
massive stars: a steeper IMF slope in the field would be 
manifested as a smaller number of field OB stars, yet the 
slope of the distribution tends to be steeper than ex- 
pected, rather than flatter. If the field star IMF is indeed 
steeper than in clusters, then the clustering law must also 
steepen substantially for the smallest clusters, in such a 
way as to compensate for a steep IMF in our data. Fur- 
ther investigation is necessary to resolve this issue. 
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Fig. 8. — Monte Carlo models for evolution in the log TV* distribution, assuming that all objects were created simultaneously. Models for 
ages of 0, 15, and 25 Myr are shown, with fitted power-law slopes. As in Figure 5, the linear power-law exponents are shown, while the fitted 
slope is —/3 + 1 in logarithmic space. 
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Fig. 9. — Spatial distribution of isolated field stars for the OB sample (plus symbols), and O-star sample (cross symbols). Scale and axes 
are as in Figure 3. 



In the meantime, the data appear to be broadly consis- 
tent with the simpler scenario of, simultaneously, a univer- 
sal IMF and universal clustering law given by equation 1. 
As emphasized by McKee & Williams (1997), the total OB 
star population, and thus the fraction of isolated field mas- 
sive stars, is driven by iV* )Up with an inverse logarithmic 
dependence. For maximum iV* )Up ranging between 10 and 
10 6 , equation 10 yields a fraction of field OB stars ranging 
from 35% to 7%, respectively. 

The clustering law in equation 1 has important conse- 
quences for feedback, implying that the interstellar poros- 
ity caused by the formation of superbubbles and supernova 
remnants has equal relative contributions from objects of 
all sizes (Oey & Clarke 1997). And, since the Stromgren 
volume V s oc TV*, it also implies that a strong majority 
of the field massive stars likely contribute to ionizing the 
diffuse, warm ionized medium, which constitutes about 
40% of the total Ha luminosity in star-forming galaxies 
(e.g., Walterbos 1998). Our result is quantitatively con- 
sistent with the result of Hoopes & Walterbos (2000) that 
field OB stars can power 40% ± 12% of the warm ionized 
medium in M33, where the fraction of field OB stars is 
likely around 15%. 

6. CONCLUSION 

We find no evidence that the field massive stars in the 
SMC are formed by a fundamentally different star-forming 



process. Rather, we find that the continuous power-law 
distribution in iV» down to N t — 1 strongly suggests that 
the star-forming process is continuous from rich clusters 
to poor groups, apparently for all ensembles that form OB 
stars. The data are consistent with the model that single, 
field OB stars are usually the most massive member of a 
group of smaller stars, as expected from the universal 
distribution (equation 1). These results are consistent with 
the simultaneous existence of a universal IMF and univer- 
sal N~ 2 clustering law. These joint universal power laws 
imply that field OB stars constitute roughly 35% to 7% of 
the total massive star population, with an inverse logarith- 
mic dependence on A r „ jUp of the the most populous cluster. 
Thus, the fraction is dependent on galaxy size and/or star 
formation rate. The contribution of these field stars to 
the ionized volume in the warm ionized medium is likely 
to scale according to their relative fraction. The universal 
clustering law also implies equal relative contributions by 
superbubbles of all sizes to the interstellar porosity. 
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Walterbos, Cristiano Porciani, Todd Small, and especially, 
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